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Abstract 

Pairing correlations in symmetric nuclear matter are studied within a rela- 
tivistic mean-field approximation based on a field theory of nucleons coupled 
to neutral ( a and w ) and to charged ( q ) mesons. The Hartree~Fock and the 
pairing fields are calculated in a self-consistent way. The energy gap is the 
result of a strong cancellation between the scalar and vector components of 
the pairing field. We find that the pair amplitude vanishes beyond a certain 
value of momentum of the paired nucleons. This fact determines an effective 
cutoff in the gap equation. The value of this cutoff gives an energy gap in 

agreement with the estimates of non relativistic calculations. 
PACS number(s): 21.65.+f, 21.30.Fe, 21.60.Jz 
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I. INTRODUCTION 



In the investigation of superfluidity or superconductivity phenomena in systems of 
strongly interacting fermions a crucial role is played by the gap equation of the BCS theory 
||l||. For infinite homogeneous systems the kernel of the gap equation is usually a slowly 
decreasing function of the quasiparticle momentum. Even if the sum over momenta con- 
verges, contributions from terms involving large momentum components of the interaction 
and quasiparticle energies very far from the Fermi surface, are not negligible [0, despite 
the fact that the gap energy may result to be a small fraction of the Fermi energy. Con- 
trary to what happens for electrons in metals, for nuclear or neutron matter, there is no 
natural cutoff over momentum. Consequently for nuclear matter a relativistic treatment is 
desireable. In the present work we investigate pairing correlations in cold symmetric nuclear 
matter by taking into account relativistic effects in a consistent way. A first study in this 
direction, for ^5*0 pairing in nuclear matter, has been already performed by H. Kucharek and 
P. Ring within the framework of the Quantum Hadrodynamics theory of Ref. ( QHD ). 
These authors, though introducing some simplifying assumptions ( e.g exchange terms of 
the interaction have been neglected ), have shown that pairing correlations and partcle-hole 
correlations can be described on the same foot within a relativistic treatment. In our opinion 
the consequences of a relativistic quantum approach have not been thoroughly exploited in 
Ref. IP]. For instance, the amplitudes of the nucleon field in the quasiparticle states have 
been assumed to be proportional simply to the four-spinors of the Hartree approximation. 

In a recent paper a more complete relativistic treatment of the Hartree-Fock- 
Bogoliubov approximation to QHD has been presented. The authors of Ref. ^ have shown 
that the pairing field has large scalar and vector components. The actual value of the gap 
in the excitation spectrum is the result of a strong mutual cancellation between these com- 
ponents of the pairing field. In the present work we obtain similar results, in particular the 
expression for the quasiparticle energy practically coincides with that of Ref. 0. However, 
though the starting point of our treatment is the same as in Ref. i.e. a relativistic gen- 
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eralization of the Gorkov scheme 0, we do not need to introduce any particular "ansatz" 
for an effective single-particle Lagrangian, as done in Ref. 0. Moreover, we show that a 
more consistent treatment, from a relativistic point of view, brings out a new interesting 
feature: the pair amplitude vanishes beyond a certain momentum of the paired nucleons. 
The value of this momentum depends on the nuclear density and on the strength of the 
self-consistent pairing field. As a consequence in the relativistic gap equation a natural 
cutoff over momentum occurs. In the previous works of Refs. it has been necessary to 
introduce such a cutoff in an arbitrary way, in order to obtain a satisfactory agreement with 
the current estimates of nonrelativistic calculations p|j9|-[T^ for the energy gap. 

We base our approach on the version I of QHD ( QHD-I ) ||^ , we include also the charged 
Q meson field in a phenomenological way similar to the MFT approach to QHD-II of 
Ref. . For simplicity we neglect the pion field. Adding the pions would produce only small 
quantitative differences in our results p. 

Here we are mainly interested in studying new effects introduced by a relativistic ap- 
proach to superfluidity. First we perform our calculations by using the same approximations 
already introduced in Refs. [0] and ||^ for studying collective modes and response functions 
of nuclear matter. These approximations amount to neglecting finite-range and retardation 
effects in the exchange of mesons between nucleons. Within this approach we obtain simple 
expressions for the relevant quantities, where we can easily appreciate the role played by the 
various ingredients of the theory. In spite of the approximations introduced, the essential 
features of a relativistically covariant treatment are retained. 

The value of the energy gap obtained with the approach just mentioned is much larger 
than the value predicted by nonrelativistic treatments P,|9|-[l4|. A better quantitative agree- 
ment can be achieved by taking into account the finite range of the nucleon-nucleon interac- 
tion. Retardation effects instead do not play a significant role. We evaluate the finite-range 
effects by using an iterative procedure and show that the first iteration gives already a suf- 
ficiently accurate approximation. 
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II. FORMALISM 



In the model adopted here the nucleons are coupled to neutral mesons ( a and uj ) and 
to the charged vector meson, q. According to the usual procedure employed in studying 
superconductivity we consider ensembles with indefinite number of particles. Therefore we 
add to the Lagrangian the term /i ?/'(x)7°-?/'(a;) , where /i is the chemical potential ( t\}{x) 
denotes the 8-component nucleon field ). In the end the value of the chemical potential 
will be determined by fixing the average baryon density. With these ingredients the field 
equations are |^ 

{idxl^ - gvVx[x)-i^ - (7,Ba(x) ■ r7^ + ^7° - M + gs^[x))^[x) = , (2.1a) 
(dxd^ + ms^Mx) = gsH^)H^) , (2.1b) 
OxW^^ix) + my V(x) = gvip{x)YHx) , (2.1c) 

9aL^"(x) + m/B^(x) = g,ilj{x)YT^{x) , (2.1d) 

where W^^^x) = d^V^x) - d^V^{x) and L^^(x) = d^B''{x) - a^B^(a;) . The quantities 
$(x), V^lx) and B'^(x) represent the scalar, vector and charged vector fields, respectively. 
In the present paper we are concerned with a uniform system at equilibrium, it is convenient 
to study such a system through the one-particle density matrix in four-momentum space: 

= J^Jd'R e-'^-^ < 0| : ^,(|)^,(-|) : |0 > , (2.2) 

where a and (3 are double indices for spin and isospin. The dots denote normal ordering 
and |0 > is the correlated ground state. Since we consider only symmetric nuclear matter 
the matrix F{p) is diagonal and degenerate with respect to isospin indices. 

By using the procedures outlined in Ref. [0, the following equation for F{p) can be 
derived from Eqs. ( p.l| ) 
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~ ^^(2^ / < °' ^ )[7V]„.^,(-|) ■ B; 

+ 9sj^, J d'R e-'^-^ < 0| : ^.(f ) : |0 >= . (2.3) 

Following Refs. and we neglect derivative terms in Eqs. ( p.lb| ,c,d) so that the 
meson field operators are simply given by 

$(x) = ■^iIj{x)iIj{x) , 
V\x) = ^^(xh^'^ix) , 

B\x) = ^iPixh^T^ix) . (2.4) 

This approximation simplifies calculations considerably, however it neglects retardation and 
finite-range effects in the exchange of mesons between nucleons. Nevertheless, because of the 
small Compton wavelengths of the heavy mesons with respect to the internucleon spacing 
in ordinary nuclear matter, the approximation ( |2.4| ) appears to be reasonable. Clearly, for 
pions this approximation would not be justified. 

After substituting in Eq. ( |2.3D the expressions ( ^.4] ) for the meson field operators we 
obtain an equation which contains expectation values of products of four nucleon field op- 
erators: 

< 0| : M^)M-§)M-§)M-§) : |0 > . (2.5) 

Following Gorkov these quantities are approximated by two-fold products of expectation 
values: 



< o| : M§)M~)n~)M~) ■■ |o > 

< 0| : M^)M~) ■■ |0 X 0| : ^^{-^)M~) : |0 > 

. , R^ ,„ „, T / R^ , R^ 



- < 0| : ^^(-)V;,(--) : |0 X 0| : ^,(--)^,(--) : |0 > 
+ < 0| : ^,(^)V;,(-^) : |0 X 0| : V;,(-^)^.(-^) : |0 > . (2.6) 



The first two terms correspond to the Hartree-Fock field. Beside these terms, Eq. ( p.6|) 
contains the product of elements of the anomalous density matrix. A nonvanishing value of 
this product implies the presence of the superfiuid phase. 
Introducing the pair amplitude 

DMp) = ^Jd'R e-'^-^ < 0| : : |0 > , (2.7) 

and the pairing field 

= J d^P D^p{p) =< 0| : Mx)Mx) : |0 > , (2.8) 
we obtain the following equation for the one-particle density matrix 

( Pa7^ - fvPBl' + /^7°)^(P) - (M - fsPsMp) 

- fvl^l^lD{p) - /,7VA ■ T^ilD{p) + fsl^D{p) = , (2.9) 
where the effective coupling constants, fs and /y, are given by the combinations 

fs = lfs + lfv + lf,, fv = lfs + lfv + lfe, (2.10) 

with fs = ( gs/ms )^ fv = ( gv/my ^ and = ( gg/m^ y. The matrix A in (^is 
conjugate to the pairing field: 

A,/3 = (7°At70)^^ = < 0| : M^)Mx) : |0 > ■ (2.11) 

The scalar density ps and the baryon density pb are given by 

ps=2Tr J d^pF{p) , 

PB=2Tr J d^p-f'^Fip) , (2.12) 

with the factor 2 coming from isospin degeneracy; the traces are taken only over the spin 
states. 

By repeating the same procedure that leads to Eq. (|2.9D, we obtain for the pair amplitude 
D{p) an equation coupled to Eq. (|2.9|) , that reads 
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( + /yPB7° - /i7°)^(p) - (M - ~fsPs)D{p) 

- fvj'A^lFip) - /,7VA • t^^IHp) + fs^Hp) = , (2.13) 
In Eqs. ( p.9| ) and ( p.l3[ ) the exchange contributions to the mean field have been taken 



into account through the definition of the effective couphng constants of Eq. ( p.lO| ). 

Before turning our attention to Eqs. ( p.9|) and (|2.13|) , we discuss the symmetry and tensor 
properties of the pairing field. Here we consider only ^Sq pairing of nucleons in symmetric 
nuclear matter. In this case the pairing field A is symmetric and degenerate with respect 
the isospin indices. Therefore we can consider only the isoscalar component of A and hence 
omit the isospin variables. Now A is a 4 x 4 matrix in spin space like the matrices F{p) and 
D{p). In spin space the matrix A is antisymmetric, hence it can be decomposed as 

A = Asa^^' + Aps(7°2 + Ao7^7' + Ai7^ + A27S° + A37' , (2.14) 

the subscripts are chosen according to the tensor properties of the various terms. The 
transformations properties of A can be derived from the basic transformation law of the field 
operator %l){x). In detail, under infinitesimal Lorentz transformations ( A'*''^ = g'^^ + e^'^ ) the 
field A transforms according to 

A' = A + '-e^'^iKaxu + ^ta/A) , 

while for space inversion 

A' = 7°A7° . 

From these equations we can see that the four quantities (Aq, — iAi, A2, iAs) represent a 
four-vector and the remaining components A^ and Aps are scalar and pseudoscalar quan- 
tities, respectively. Since we are considering a homogeneous system at rest, only the com- 
ponents As and Aq in Eq. (|2.14| ) can differ from zero: 



A = A5(t1^ + A07V • (2.15) 

Moreover, because of the invariance of the equilibrium state under time-reversal the com- 
ponents As and Aq can be assumed to be real and the matrix A of Eq. ( p.ll|) becomes 



A = A . (2.16) 

From Eqs. ( p.9| ) and ( p.l3| ), after some algebraic manipulations, we obtain two separate 
equations for the density matrix F(p) and for the pair amplitude D(p): 

{R^{ph^-M{p))F{p) = (2.17) 

and 

(R^i-pH - M{~P)) Dip) = . (2.18) 
The components of the four-vector R^ip) are 
R= p , 

Roip) = + /i + ^ [g'Alfi - fAlpo + fgAsAoM*) , (2.19) 

and the mass term M{p) is given by 

M(p) = M*--^{fAlM*-fgAoAs{po-fi)) , (2.20) 

Here the coupling constants g and / are given by the combinations 

9= fs-^{fv + f-0) , f= fs + 2{fv + f,) , (2.21) 

and the quantity /C(p) is expressed as 

/C(p) = {Po - fif -El + fAl- g' Al , (2.22) 

where Ep = (p^ + M*^)^/^, while fi = ^ — fvpB is the effective chemical potential and 
M* = M — fsps is the effective nucleon mass. 

Equation ( p.l7| ) tells us that the matrix F{p) can be put in the form 



F{p)= F{p)+rF^{p) , (2.23a) 

with 



and F{p) a scalar quantity. 

Equations ( |2.17| ) and ( p.l8| ) contain the components of the pairing field as parame- 



ters, these components must be determined self-consistently. This can be done by solving 
Eq. ( |2.13| ) with respect to D(p) with the aid of Eqs. ( |2.23| ). For the two nonvanishing 
components of the pair amplitude 



D{p) = Ds{p)a'' + Doiph'^^ 



we obtain 



Dsip) = ^ (g^sfiFoip) - fAoipoF{p) - M*Foip))) , (2.24a) 
IC{p) 

Doip) = ^ {g^smv) - fMPoFoip) - M*F{p))) . (2.24b) 

Now we derive the energy spectrum using our approach. Substituting in Eq. ( |2.17] ) the 
formal solution given by Eqs. ( |2.23| ), we can see that the components of the four-vector 
Rjiip) must satisfy the constraint 

R^R^ = m\p) . 
This equation, for a fixed p, determine the allowed values of po 

2 ifi'El + g'fAlAl + fAy - 2fg A^Ao fiM* f' , (2.25) 

which correspond to the energies of elementary excitations referred to the effective chemical 
potential /i. The upper sign refers to elementary excitations of the Dirac sea. Actually in 
the absence of the pairing field the energies of these excitations become 

Po = {Ep + /i) . 
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We neglect contributions from the Dirac sea and consider only the energy value given by 
the lower sign in Eq. ( |2.25| ). This defines the energy Sp of a quasiparticle with momentum 
p in the superfiuid phase. 

We turn now to the explicit evaluation of the matrices F(p) and D{p). By inserting a 
complete set of energy eigenvectors in Eq. (|2.2|) we have 



^-P(P) = J^ ld'Re'^-^Y.^(Po + En)<0\M^)\r^><ri\M-^)\0> , (2.26) 

where \n > represents a quasiparticle state \n >= |£^p„,Pn, A„ > with spin label A„. This 
equation, toghether with Eq. (|2.17|) , shows that the one-particle density matrix F{p) can 
be put in the form 

Fa,{p) = j^sT.9l{p)4'\pWx\p)9xiPmPo + Sp) , (2.27) 

where the spinor ux{p) obeys the equation 

(i?^(p, -£p)Y - M{-£p)>xiP) = . (2.28) 

With the normalization Ux{p)u\{p) = 1, the quantity gx{p) represents the probability am- 
plitude of finding a hole with momentum — p and spin label —A in the quasiparticle state 
\£p,P, X >. 

For the pair amplitude D(p), with the same procedure and with the aid of Eq. ( p. 181) , 
we obtain the expression 

Dapip) = j^,T.i-^y9Up)uTipWx\p)UPmPo + Sp) . (2.29) 

The equation for the spinor v\{p) can be derived from Eq. ( p.l8| ) and reads 

{R,{-p, 8p)r - m^))Mp) = . (2.30) 



We choose for vx{p) the same normalization as ma(p)- In Eq- ( p.30| ) the quantity gx{p) is 



the probability amplitude of finding a particle with momentum p and spin label A in the 
quasiparticle state |£^p,p, A >. With the choice of phase made in Eq. (p.29|) the product 
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(7* a(p)^a(p) is real and independent of A, as can be seen by substituting in Eq. ( |2.13|) the 
expression ( |2.27|) for F{p). 

It is implicit in our approach that a quasiparticle state is a superposition of one-particle 
and one-hole states, so that the following normalization condition holds 



I^7a(p)P + |^a(p)|'= 1 . 



(2.31) 



Moreover, since we are considering a homogeneous and isotropic system, |5'a(p)| and |^a(p)| 
do not depend on A. 

The specific form of the spinors ux{p) and vx{p) is determined by the signs of Ro{p) 
and M{p). Explicit calculations show that -Ro(=t^^p) remains positive for any value of |p|, 
whereas the mass term M{—Sp) of Eq. ( |2.28| ) has a peculiar behaviour. It is a monotonically 
decreasing function of |p|, it takes negative values for |p| larger than a certain value Pc and 
becomes infinitely negative at a finite value of |p| > Pc- The mass term M{Sp) of Eq. (p.30|) 
instead is always positive and almost constant. 

For IpI < Pc the solutions of Eq. ( p^.28D are given by 

_ / XA \ 

, , ri?n(-0 + M(-fJii/2 
«a(p) = 



2Rq{—£p) 

^ Ro{-Sp) + M{-Sp 
whereas for |p| > pc we have to choose the spinors 



( 



«AlPj 



\M{-£p)\ 



1/2 



a p 

tt^^ttXa/ 



Xa\ 



i?o(-^p) + |M(-^p)| 



V 



Xa 



as solutions of Eq. (|2.28| ) . The solutions of Eq. (|2.3CI| ) are given by the spinors 

_ ( XX \ 

va(p) = 



V 



<rp 



Roi£p)+M{£p) 

In equations above x\ denote the usual two-component Pauli spinors. 



(2.32a) 



(2.32b) 



(2.33) 
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We remark that the energy spectrum of excitations of the Dirac sea remains well sepa- 
rated from the quasiparticle spectrum. In fact the mass terms in Eqs. ( |2.17] ) and ( p.l8|) for 
the antiparticle case, are always positive for any value of |p|. 

By substituting the spinors ( p.32b| ) and ( p.33| ) in Eq. ( p.29| ), we can see that for |p| > pc, 
the two components 

and 

Do(p) = ^Tr7V^(p) 

vanish. This implies that -fo(p) also vanishes for |p| > pc, see Eqs. ( |2.24j) . The occupation 
number of particles in the correlated ground state displays a discontinuity. 

From Eqs (|2.24|) with the aid of the normalization condition (p.31|) , we can determine the 
amplitudes g\{p) and ^a(p), and the matrices F{p) and D(p), as functions of the parameters 
M*, /i, As and Aq. These parameters can be calculated by solving the four coupled equations 
that are obtained by fixing the baryon density 



2 1 

3 ^ 



pl=^j rfpFo(p) , (2.34a) 



and using the self-consistency relations for the effective nucleon mass 



M*= M-8fs I dp ) / Fo(p) (2.34b) 



and for the components of the pairing field 



^s=2 f ^ , (g^sfiFoip) - /Ao(-^pF(p) - M*Fo(p))) , (2.34c) 

^0=2/ r/^ (9\s^F{p) - /Ao(-^pFo(p) - M*F{p))) . (2.34d) 
The momentum pc, for which -Fo(p) vanishes, plays the role of an effective cutoff in 



the integrals (|2.34|) . This fact avoids introducing an artificial cutoff to make the integrals 
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( p.34c| ) and (|2.34d| ) converge. The dependence of pc on As and Aq amounts to a further 



self-consistency relation. 

Equations ( |2.34c| ) and ( p.34d| ) replace the gap equation of the nonrelativistic case. In 



the nonrelativistic limit F{p) -fo(p); hence these two equations become identical and the 
two components of the pairing field coincide: A5 = Aq. Moreover, the quasiparticle energy 
Sp acquires the usual expression of the BCS theory 



where A = A^ = Aq. 



III. RESULTS 



In this section we investigate the solutions of Eqs. ( 2.34|) . These equations have been 



derived by starting from the assumptions ( p.4|) , which neglect the finite range of the nucleon- 
nucleon interaction. 

Concerning the coupling constants fv and fs we choose their value so as to reproduce 
the binding energy ( 15.75 MeV ) of saturated nuclear matter with a Fermi momentum of 
1.42 fm~^ ( see Ref. ). These values are: 

fs = 2.37-10-^MeV~\ fv = lA^-W'^MeV-^ . 
For the coupling constant we have taken the value determined by the g 2tt decay 

f, = 1.55-10-^ MeV-'^ . 
Then the effective constants / and g for the pairing field are 

/ = 5.58-10~^Me\/-^ g = -Am-lO-^MeV'^ . 

The relevant quantities for the quasiparticle energy spectrum are f Aq and gAs- In 
Fig. |I| these quantities are displayed, together with their sum, as functions of the Fermi 
momentum. The sum f Aq + g As approximatively reproduces the gap in the quasiparticle 
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energy spectrum. Figure |l| shows that a superfluid solution of Eqs. ( p.34|) is present in the 



range of densities corresponding to pp < 1.35 /m^^. Moreover one can see that / Aq and 
g As separately are very large with respect to their sum, i.e. the gap in the quasiparticle 
spectrum is determined by difference between two large and not very different numbers. 

In Fig. 1^ the excitation spectra of the superfluid and normal phases are displayed as 
functions of |p| for pp = 0.9 /m^^. A numerical analysis shows that for |p| >~ 1.9 /m^^ 
the quasiparticle energies lie slightly below the excitation energies of the normal phase, 
Ep — jl. This fact gives rise to the cutoff in the integrals (|2.34| ). Since for a sufficiently high 



value of IpI the quantity KL{p) of Eq. (|2.22|) may vanish, the second term of M(p), Eq. ( p.20D , 



which is positive, can become larger than M*. The cutoff pc is given by the value of |p| for 
which the r.h.s of Eq. ( |2.20|) vanishes. 



In Fig. we show both the sum / Aq + A5 and the energy gap as a function of the 
Fermi momentum. We can see that, though remaining small, the difference between these 
two quantities increases with pp. This is because relativistic effects become more important 
with increasing density. In the non relativistic limit the energy gap and the sum / Aq + (? A5 
coincide. Then, the difference between these two quantities could give an insight about the 
relevance of relativistic effects. In the region around pp = 0.9 fm^^ where the energy gap 
takes its maximum value, this difference is not very large, only about 10%. However we 
remark that the occurrence of a cutoff in the integrals for the pairing field is obtained only 
using a relativistic expression ( Eq. ( p.25|) ) for the quasiparticle energy. 

Figure ^ shows that the pair amplitude is rather sharply peaked about the Fermi mo- 
mentum. This important feature is the basis of the approximation, that we use in the next 
section to take into account the finite range of the nucleon-nucleon interaction. 

For the two values of the Fermi momentum for which superfluidity disappears and the 
energy gap becomes maximum, our results qualitatively agree with previous treatments, both 
nonrelativistic PJ9|-p!^ and relativistic Instead for the most relevant quantity of the 

superfluid phase, the energy gap, our calculations yield values that are twice the generally 
accepted estimates. The value of the sum f + g, which plays the role of an effective coupling 
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constant for interacting paired nucleons is too large. It is worth noticing that the role oi f+g 
in determining the energy gap is enhanced by a cooperative effect due to the self-cosistency 
constraint for the cutoff pc- if the gap becomes larger the value of pc increases, then the 
contribution to the pairing field of the integrals ( |2.34c| ) and ( |2 .34011 ) also increases, giving 
rise to a larger gap. 

IV. FINITE RANGE EFFECTS 

In this section we evaluate effects due to the finite range of the nucleon-nucleon inter- 
action. For the effects of retardation in the propagation of the meson fields we will give 
only an estimate of their magnitude and argue that these effects do not play an important 
role. For simplicity we take into account the finite range of the effective interaction between 
paired nucleons only when calculating quantities that are relevant to the superfiuid phase. 
In deriving the Hartree-Fock field, instead, we retain the approximation of Eqs. ( p.4| ). Thus 
for quantities containing the pairing field, instead of the approximated expressions ( p.4| ) we 
introduce the formal solutions of Eqs. ( p.lb| ,c,d): 

"^{x) = 9s J d'^yV{x - y)i){y)^{y) , 

Va(x) = gv j - y)i^{y)i'i^{y) , 

BA(a;) =g,j d'y vi'^ix - ymyh^r^iy) , (4.1) 

where V{x — y), V^^J{x — y) and V^x/ii^ " y) the propagators of the meson fields. 

We explicitly derive the equations for the contributions of the scalar field alone, for the 
vector fields only the final results are reported. 

In Eq. ( |2.6|) the term containing elements of the anomalous density matrix acquires the 
form: 

j d'yViy - f ) < 0| : M^Wv) ■ |0 >< 0| : : |0 > , 

then in Eq. ( |2.9| ) the term fsAD(p) is replaced by 
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gl I d'qV{p-q)fD\q)^'D{p) . (4.2) 
Here V{p — q) is the Fourier transform of the propagator P(a; — y). 



Analogously, in Eq. ( p.l3|) we have to make the substitution 



/5 A F{p) ^ glj d'q V{p - q)D{q) F{p) . (4.3) 

Assuming that retardation effects are negligible, we can put Po — qo = in evaluating 
the integrals of Eqs. ( |4.2|) and ( [4.3|) . Thus these integrals depend only on p. This fact 
greatly simplifies calculations. In particular, the quasiparticle energies can still be expressed 
in closed form; it is sufficient to replace in Eq. (|2.25| ) the terms gAs and /Aq with the 
analogous p-dependent quantities. 

We introduce now a further and more fundamental approximation. In the previous 
section we have shown that the components of the pair amplitude D{q) are strongly peaked 
about Pp. The width of this peak is much smaller than the range over which the meson 
propagators present an appreciable variation. In fact this range is typically of order ~ 
ms^mv ^ttIq. For this reason we expect that the values of the integrals (|4.2| ) and ( ^^) can 
be given by 

glV{\^lpp) I d'q^'D^q)^' = glVi\p\,PFh'^^7' , (4.4) 

gln\p\,pp) J d'qD{q) = glV{\p\,pp)A , (4.5) 

with a satisfactory approximation. Here V{\p\,pf) stands for the average over the directions 
of q. Thus the quantity /5A in Eqs. (|2.9|) and ( p.l3|) is simply replaced by 

fsasilplpp)^ , (4.6) 

where the factor 

4:Pf\p\ \m''s + {\p\ - PfY J 
represents the finite-range correction to the contribution of the scalar meson. 
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As far as the vector mesons are concerned, with some algebraic manipulations we can 
see that the matrices fy'j^A'jJ and Jv'J^tA ■ T'jJ in Eqs. ( p.9| )and ( |2.13| ) must be replaced 
by 

4/y(^ay(|p|,PF) + \)As^'' - 2fv{\ay{\p\,pF) + ^)Ao7V (4.7a) 
for the u} meson, and by a similar expression for the g meson, 

4/.(^a.(|p|,PF) + \)Asa'' - 2f/-aMPF) + ^)Ao7V • (4.7b) 

Here 

2 

«v(|p|,Pf) = ml.V{\p\,pF) = - — f- In 



aeilPl^Pp) = mfD{\p\,pF) = - — ^ In' ^ 



ml + (IpI 


+ Pf? 


mf. + (IpI 





The correction factors are given by the bracketed terms in Eqs. (f4.7|) . 

Summarizing, in the present approximation the finite range of the nucleon-nucleon inter- 
action is taken into account by replacing in all the equations of Sect. II the pairing coupling 
constants / and g, with the combinations 

a/(|p|,PF) = fsas{\vlPF) +/y(ay(|p|,PF) + 1) + ^(ae(|p|'Pi^) + 1) (4.8a) 

and 

ag(|p|,PF) = fsas{\vlPF) - /y(3av(|p|,PF) + 1) - /e(3ag(|p|,PF) + 1) (4.8b) 
respectively. 

The pairing coupling constants, separately considered, are not much affected by the 
finite-range corrections. In fact, the values of a/(|p|, pi?) and ag(|p|,pir) for |p| = p/, with 
0.6 /m~^ < Pf < 1.0 /m^^, are smaller than / and g by about 10% and 2%, respectively. 
However, for the sum / + (?, that practically determines the magnitude of the energy gap, the 
correction is more important. In the same range of pf the value of afipF.pp) + CLgipF^PF) 
is about |(/ + (7). 

17 



The most important effects of the finite-range corrections are a substantial quenching 
of the pairing field and a reduction of the energy gap by an overall factor three or more. 
This is shown in Figs. |^ and ^ where the calculated components of the pairing field and 
of the energy gap are shown as functions of pi;'. Moreover, with respect to the zero-range 
approximation the domain oi pp where the superfiuid phase can arise, is narrowed and the 
maximum of the energy gap is shifted towards a lower value of pp, 

In Fig. the quasiparticle energy together with the excitation spectrum of the normal 
phase is displayed for pp = 0.8/m~^, where the energy gap takes now its maximum value. 
The value of |p| where the two curves cross, corresponds to the cutoff pc- A numerical 
analysis shows that Pc ~ 1.7 fm~^. 

Finally, in Fig. | the components of the pair amplitude are shown as functions of |p| for 
Pf = The finite-range corrections makes the pair amplitude even more peaked 

about Pf- This fact further justifies the approximation expressed by Eqs. ( [4.4|) and ([4.51). 
That approximation simplifies calculations substantially. In fact, the basic quantities, i.e. 
the components of the pairing field, that are determined self-consistently, are still two con- 
stants. We have assessed the reliability of this approximation by successive iterations of 
Eqs. ( |4.2| ) and ( [4. 3) ), starting from Eqs. ( |4.4] ) and ( [4.5| ). The correction coming from the 
first iteration amounts to a few percent, while the second iteration does not substantially 
modify the first-order results. 

We have evaluated also the order of magnitude of retardation effects, replacing £p and 

instead of po and go in Eqs. ( [4.2|) and ([4.31). The change is less than one percent. 



V. SUMMARY 

We have investigated in the framework of a relativistic model, the possibility for the onset 
of a superfiuid phase in symmetric nuclear matter and the features of this phase. We have 
derived equations for the relevant quantities of the superfiuid phase by making a relativistic 
generalization of the scheme introduced by Gorkov to study superconductivity in electron 
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systems. In this scheme nuclear matter is described as an ensemble of quasiparticles moving 
in the Hartree-Fock field plus a self-consistent paring field. At first both the Hartree-Fock 
field and the pairing field have been treated on the same foot by using an approximation 
where the finite range and the retardation of the meson propagation between nucleons have 
been neglected. Then, we have improved our approach by introducing finite-range effects 
for the quantities pertaining to the superfluid phase. 

In the relativistic treatment of the pairing process for a system at rest the pairing field 
has two different components: a Lorentz scalar and the time-component Aq of a four- 
vector. The different behaviour of the two components under Lorentz transformations must 
be properly taken into account in the relativistic hydrodynamics of nuclear systems in the 
superfluid phase. 

In the non relativistic limit A^ and Aq coincide. In our approach they are slightly dif- 
ferent. This does not mean, however, that relativistic effects are negligible. Actually, the 
expression for the quasiparticle energy derived in our calculations, differs from the nonrel- 
ativistic expression of the BCS theory in an essential way. The quasiparticle energy given 
by Eq. ( |2.25| ) or by its analogous expression that includes flnite-range corrections, displays 
the salient feature that it can be less than {Ep — jl) beyond a certain value of |p|. This 
determines the occurence of a cutoff in the relativistic gap equations ( p.34c|) and (|2.34dD , 
that does not appear in the analogous nonrelativistic expression. This cutoff removes the 
contributions of high |p| components of the interaction to the gap equation, thus allowing 
the use of nucleon-nucleon interactions that are only slowly decreasing for high values of 
|p|. This fact reduces the energy gap appreciably. 
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FIGURES 

FIG. 1. Components of the pairing field ( times the respective couphng constants ) as a function 
of pF- The sohd hne and the dashed Une correspond to the vector and the scalar components, 
respectively. The dotted line gives the sum of these two quantities. 

FIG. 2. Quasiparticle energy ( solid line ) and single-particle energy Ep — Jl ( dashed line ), 
iov pp = 0.9fm~^. 

FIG. 3. Energy gap ( solid line ) together with the sum /Aq + gAs ( dashed line ) versus pp- 

FIG. 4. Components of the pair amplitude ( times (27r^)^ ) as a function of |p| for 
Pf = 0.9 fm~^. The solid and the dashed lines correspond to the vector and scalar components, 
respectively. 

FIG. 5. Components of the pairing field at the Fermi surface as functions pp- Finite-range 
corrections are included. The vector component ( solid line ) and the scalar component ( dashed 
line ) are respectively multiplied by and ag ( see Eqs. ( |4.8| ) ). 

FIG. 6. Energy gap ( solid line ) versus pp together with the sum o/Aq + a^A^, evaluated at 
the Fermi surface ( dashed line ). The finite-range corrections are included. 

FIG. 7. Same as Fig. 2 with finite-range corrections and for pp = 0.8 fm~^. 

FIG. 8. Same as Fig. 4 with finite-range corrections and for pp = 0.8/m~^. The two compo- 
nents of the pair amplitude are practically indistinguishable at this density. 
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Fig. 4 




Fig. 8 



